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Simulation  of  the  Femtosecond  Optical  Response  of  a  Solute  in  Water 

Laurence  E.  Fried.  Nathan  Bernstein,  and  Shaul  Mukamel 
Department  of  Chemistry.  University  of  Rochester.  Rochester.  Sew  York  14627 
(Received  Q  August  199!) 

We  present  microscopic  semiclassical  simulations  of  the  femtosecond  nonlinear  optical  polarisation 
iPlx>)  of  a  solute  in  water.  P{"  is  dominated  at  short  times  by  mtermolecular  vibrations  that  cannot  be 
accounted  for  by  the  Bloch  equations.  The  applicability  of  the  Bloch,  stochastic,  and  Brownian  oscillator 
models  to  the  interpretation  of  photon  echo  and  pump-probe  measurements  is  discussed.  We  find  that 
only  microscopic  semiclassical  simulation  and  the  Brownian  oscillator  model  can  account  for  the  essen¬ 
tial  features  of 

PACS  numbers:  42.65.— k.  78.20  — e 


Femtosecond  nonlinear  optical  spectroscopy  of  dyes  in 
liquids  cannot  be  interpreted  using  the  Bloch  equations. 
The  system  has  a  multitude  of  relaxation  time  scales  and 
in  general  they  cannot  be  considered  to  be  either  infinitely 
fast  (homogeneous  broadening)  or  infinitely  slow  (inho¬ 
mogeneous  broadening)  on  the  experimental  time  scale. 
This  state  of  affairs  was  recently  highlighted  by  photon 
echo  (PE)  measurements  in  polar  liquids  1 1.2].  The 
nonexponential  decay  of  the  PE  found  in  these  experi¬ 
ments  indicates  that  the  Bloch  level  of  description  is 
inadequate.  The  origin  of  the  non-Bloch  behavior  in 
liquids  is  the  continuous  distribution  of  solvent  time 
scales,  rather  than  the  strong  field  effects  related  to  the 
Rabi  frequency  observed  previously  I3l. 

In  this  Letter  we  present  a  microscopic  semiclassical 
simulation  {4,5}  of  the  nonlinear  polarization  P<3)  of  a 
dye  in  water.  We  reproduce  many  essential  features  of 
the  recent  experiments,  including  the  non-Bloch  behavior 
of  the  photon  echo  signal.  We  use  the  results  of  our 
simulation  to  test  the  Brownian  oscillator  model  16],  the 
stochastic  model  17,8],  and  the  Bloch  equations.  This  al¬ 
lows  us  to  determine  the  level  of  theoretical  modeling 
sufficient  to  interpret  these  recent  experiments.  Previous¬ 
ly,  photon  echo  signals  in  polar  liquids  have  been  inter¬ 
preted  using  the  stochastic  model  in  the  limit  of  an  over¬ 
damped  solvent  (1,2).  Our  simulations  show  that  inter- 
molecular  vibrations  are  very  important  in  the  short-time 
optical  polarization  of  a  solute  in  a  polar  liquid  We  find 
that  only  the  Brownian  oscillator  model  and  microscopic 
simulation  can  reproduce  the  essential  features  of  the 
nonlinear  optical  polarization  of  a  solute  in  a  polar  liquid. 

The  present  microscopic  semiclassical  approach  is  de¬ 
rived  by  taking  the  classical  limit  of  exact  quantum  ex¬ 
pressions  for  P<3)  (4,5],  It  takes  into  account  a  distribu¬ 
tion  of  solvent  time  scales  and  the  difference  between  dy¬ 
namics  in  the  ground  and  excited  states  of  the  chromo- 
phore.  The  Brownian  oscillator  model  is  a  limiting  case 
of  the  microscopic  semiclassical  theory  applicable  when 
the  ground-  and  excited-state  dynamics  are  similar.  The 
stochastic  model,  in  turn,  can  be  derived  from  the 
Brownian  oscillator  model  by  assuming  an  overdamped 
solvent  and  neglecting  the  effect  of  the  solute  on  the  sol¬ 


vent.  This  means  that  solvation  dynamics  is  completely 
missed.  Finally,  the  Bloch  equations  may  be  derived 
from  the  stochastic  model  by  assuming  that  the  solvent 
has  only  infinitely  fast  and  infinitely  slow  time  scales. 
The  four  methods  thus  constitute  a  hierarchy  of  approxi¬ 
mations. 

In  our  simulation  the  solute  is  assumed  to  have  two 
relevant  electronic  states,  |g)  (the  ground  state)  and  !e> 
(the  excited  state).  The  Hamiltonian  for  this  system  in 
the  presence  of  an  electric  field  E  U )  is  *"  W  ~VE(t). 
where 


H-\g)HK(Q)(g\  +  \e)Hf{Q)  +  o>'M 

and  V  —  •e)p<g|  +  |g)/i<e|.  Here  Hg  and  W,  arc  the  adia¬ 
batic  Hamiltonians  of  the  solute  plus  solvent  system.  The 
set  of  all  nuclear  positions  is  represented  by  the  vector  Q. 
wfg  is  gas  phase  0-0  absorption  frequency  of  the  solute 
Below,  we  neglect  the  dependence  of  the  dipole  operator 
on  nuclear  position  (Condon  approximation),  the  internal 
vibrations  of  the  solute,  and  the  solute  motion  during  op¬ 
tical  probing. 

Hg  and  He  are  of  the  form  tf, -T+ F<w  w)  + K,(w  s>. 
where  i  **g,e,  T  is  the  kinetic  energy.  K(3V  is  the 
water-water  interaction  (taken  to  be  the  same  in  the 
ground  and  excited  electronic  state),  and  F,(W  S>  is  the 
water-solute  interaction.  Both  the  water  molecules  and 
the  solute  were  taken  to  be  rigid.  Interaction  site  models 
were  used  for  the  potential  between  molecules;  i.e.,  the 
yth  atom  of  each  molecule  was  assigned  Lennard-Jones 
parameters  o,,t,  and  a  partial  charge  qt.  This  leads  to 
the  atom-atom  potential  form 


+  4 


<T,  +  <7 J 

12 

a,  +  a, 

b ’ 

1  2r>J 

2r» 

( 

where  riy  is  the  distance  between  atoms  /  and  j.  The  sim¬ 
ple  point  chargc-enhanced  potential  [9]  provided  param¬ 
eters  appropriate  to  water.  The  solute  had  one  Lennard- 
Jones  site  (^“0)  at  the  molecular  center  of  mass  with 
e  —  70  cm'1  and  o-*4  A.  The  Lennard-Jones  site  was 
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taken  to  be  the  same  in  the  g  and  c  states.  In  addition,  terms  16}:  Pn’(r  %.t  :.i , ) -P'"(i  (,r:,r , )  +  P^Ht  u.’.D ) 
the  solute  in  the  excited  state  had  two  fixed  charge  sites  P"’  has  the  following  interpretation:  The  first  pulse 
(ey  ”0)  with  opposite  charges  of  magnitude  q,  located  sets  up  an  optical  coherence  in  the  sample  During  /  the 

at  positions  (0.0,  ±0  25  A).  This  model  represents  a  typ-  evolution  (and  dephasing)  of  this  coherence  occurs  The 

ical  nonpolar  dye  molecule  with  a  large  dipole  moment  in  second  pulse  converts  the  optical  coherence  into  a  popula¬ 
te  excited  electronic  state.  tion  on  either  the  ground  (P*1')  or  excited  (P.,u)  elec- 

The  signal  S.ipiAf  i./;)  in  the  three-pulse  photon  echo  tronic  slate  During  i:  thermal  fluctuations  occur  and  the 

(3PE)  is  given  by  „S\ii*t (t  i.»;)  “/“dr  \\P *■"(/ i.r  :,r  i  )| \  system  relaxes  toward  equilibrium  This  relaxation  leads 

where  t\  is  the  delay  between  the  first  two  pulses,  i:  is  the  to  the  Stokes  shift  in  P' ' 1  The  third  pulse  converts  the 

delay  between  the  second  and  third  pulses,  and  ft  is  the  population  back  into  a  coherence,  which  dephases  during 

interval  between  the  third  pulse  and  detection.  The  ordi-  the  i\  evolution. 

nary  photon  echo  is  a  special  case  of  the  3PE:  ,S'i»|  (f|)  The  microscopic  semiclassical  approximation  for  Piu 
=*5ipt(/i,<:  “0).  P0>  can  be  expressed  as  a  sum  of  two  [4,5]  (for  pulses  that  are  short  compared  to  pure  dephas- 
- ing  and  molecular  time  scales)  is 

Pt,,,,(rJ,f:Ti)-J'  dQdPexp  (/  £  \  (<?[r'./,])df')exp  ( ~i  J\' (Q[i'\  )dt']pt(Q.P)  . 

(2) 

P<"UU2.lt)*’fdQJPexp{ifi':  +  "HQlt']  )di'}exp[-i  f\(QUVdr}Pg(Q.P) 


Here  p=E(<w-s,-p( ,(w  s).  is  the  classical  Boltzmann 
distribution  of  the  g  potential  surface.  Q\t\  is  defined  as 
the  coordinate  of  initial  condition  ( Q.P )  propagated  for  a 
period  l  with  classical  mechanics  on  the  ground-s  ate  (g) 
surface.  is  defined  as  the  coordinate  of  (Q.P) 

propagated  for  a  time  f|  on  the  g  surface,  and  then  time 
ti  on  the  e  surface. 

We  evaluated  Eq.  (2)  by  performing  a  constant-energy 
molecular-dynamics  simulation  with  100  water  molecules 
and  I  solute.  Periodic  boundary  conditions  were  used, 
with  the  box  size  adjusted  to  give  a  density  of  I  g/cm3. 
The  reaction  field  method  [10]  was  used  to  correct  for  the 
presence  of  the  long-range  Coulomb  forces.  Because  of 
the  ergodicity  of  the  fluid,  we  were  able  to  generate  initial 
conditions  (Q.P)  by  taking  samples  at  intervals  along  a 
single  trajectory  on  the  ground  potential  surface.  This 
“base"  trajectory  was  generated  by  starting  from  a  ran¬ 
dom  initial  condition  and  then  performing  a  Monte  Carlo 
equilibration  for  30000  cycles.  The  equilibrated  config¬ 
uration  was  used  for  the  initial  coordinate  of  the  base  tra¬ 
jectory.  Initial  velocities  for  the  base  trajectory  were 
chosen  from  the  Boltzmann  distribution.  The  base  trajec¬ 
tory  was  run  for  a  time  of  120  ps  using  a  time  step  of  0.9 
fs.  The  coordinate  along  the  base  trajectory  gave  @[rj. 
Every  256th  step,  a  secondary  trajectory  was  run  on  the 
excited-state  surface  for  1024  steps.  This  served  to  gen¬ 
erate  Q[t2.t  i).  The  total  length  of  the  simulation,  includ¬ 
ing  both  base  and  secondary  'rajectories,  was  600  ps. 

We  first  consider  the  ground-state  and  excited-state 
dynamics  generated  by  the  simulation.  The  autocorrela¬ 
tion  function  of  U(Q{t\),  C(/)s(C(£>(r])L'(()[0])) 
-<C((?(0]>2)  (where  ()  denotes  thermal  expectation 
value),  is  directly  relevant  to  spectroscopic  measure¬ 
ments.  Within  the  Brownian  oscillator  model  all  spectro¬ 
scopic  observables  are  controlled  by  this  autocorrelation 
function  [63.  The  mapping  onto  the  Brownian  oscillator 


model  was  achieved  by  attributing  the  classical  simula¬ 
tion  to  the  symmetrized  correlation  function  C’(/) 
+  C(—  f),  and  using  the  fluctuation-dissipation  theorem 
to  get  the  complex  C(i )  [6]  The  stochastic  model  direct¬ 
ly  attributes  the  classical  simulation  to  C'(r)  which  is 
therefore  taken  to  be  real  [6]  C(i)  was  evaluated  along 
the  base  trajectory  described  above 

Figure  I  shows  C(t)  for  the  solute  in  water.  CO)  is 
clearly  in  the  intermediate  modulation  regime,  where 
Veto)  is  of  the  same  order  of  magnitude  as  the  decay 
rate  of  C(i).  This  is  in  accord  with  recent  experimental 
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FIG.  I.  Left  scale:  solid  line,  the  autocorrelation  of  the  ener¬ 
gy  gap  CO);  dolled  line,  til  by  six  Brownian  oscillators  Right 
scale:  (ViQll.O])). 
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results  [1.21.  C(t)  shows  a  fast  (period  of  45  fs)  vibra¬ 
tional  component  of  the  solvent  relaxation  III].  The 
time-scale  separation  into  very  slow  and  very  fast  terms 
assumed  by  the  Bloch  equations  is  not  apparent  in  the  ac¬ 
tual  CO).  In  the  Brownian  oscillator  model.  CO)  is  the 
sum  of  /V  terms:  C(/)*“C(0)5I,P,-C|(f)  l6l.  Here  P,  is 
the  normalized  weight  of  oscillator  i.  C,U )  depends  on 
the  oscillator's  frequency  to,  and  damping  y,.  Shown  in 
Fig.  I  is  a  fit  to  CO)  based  on  six  Brownian  oscillators. 
The  parameters  used  (/>,.to,,y, )  were  (0.55,135.539), 
(0.18,817,277),  (0.08.602,243),  (0.07,452.439).  (0.06. 
379,420).  and  (0.04,49,752),  where  to,  and  y,  are  given  in 
cm’1.  C(0)  “4.6  x  105  cm 

Figure  1  also  shows  (L/((?l/,0])>  found  upon  excitation 
from  the  ground  state  to  the  excited  state.  Compared 
with  the  ground  state,  the  excited-state  dynamics  shows  a 
more  prominent  slow  relaxation:  while  C(f)  has  relaxed 
to  within  5%  of  its  asymptotic  value  with  500  fs. 
((,'(()[/, 01 ))  has  relaxed  only  to  within  35%  of  its  asymp¬ 
totic  value  of  =  3400  cm  within  500  fs.  According  to 
the  Brownian  oscillator  model,  (E/(()[f.0] ))  C(t ).  The 
departure  from  this  condition  demonstrates  that  the 
Brownian  oscillator  model  does  not  give  a  quantitative 
description  of  excited-state  dynamics. 

Linear  and  nonlinear  spectroscopies  provide  a  wide 
variety  of  “windows"  onto  the  ground-  and  excited-state 
dynamics  discussed  above.  PE,  3PE.  and  pump-probe 
spectroscopy  probe  different  parts  of  P<}>  and  the  under¬ 
lying  molecular  dynamics.  We  begin  our  analysis  with 
linear  optics:  The  inset  of  Fig.  2  shows  the  ordinary  ab- 
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FIG  2.  (a)  The  two-pulse  photon  echo  signal.  Solid  line, 
calculation  of  Eq.  (2);  dotted  line,  cumulant  calculationfheavy 
dotted  line,  best  fit  of  Bloch  equations  to  the  PE  signal,  (b)  As 
in  (a),  but  the  logarithm  of  the  signal  is  shown.  Inset:  Absorp¬ 
tion  and  fluorescence. 


sorption  and  fluorescence  spectra,  calculated  by  Eq  (2) 
The  frequency  to  is  shown  relative  to  Both  spectra 
are  broad  and  featureless,  and  hence  contain  little  de¬ 
tailed  dynamical  information.  The  fluorescence  displays 
a  Stokes  shift  of  about  3400  cm  " 1  to  the  red  of  the  ab¬ 
sorption  line,  while  the  Brownian  oscillator  model  pre¬ 
dicts  a  Stokes  shift  equal  to  2200  cm  ” 1 .  The  Stokes  shift 
is  missing  in  the  Bloch  equation  and  stochastic  ap¬ 
proaches. 

In  Fig.  2(a)  we  show  the  PE  calculated  with  the  Pen> 
term  of  Eq.  (2)  only  (the  quantum  expression  (6]  for  P"> 
yields  P">(r, ,/y-O./,) -P,‘ "(/ ,,r; -O.r , )].  The  initial 
rise  time  is  due  to  the  finite  inhomogeneous  broadening  in 
this  system:  in  the  limit  of  larger  inhomogeneous 
broadening  the  induction  time  decreases.  Nibbermg, 
Wiersma.  and  Duppen  [2]  have  observed  a  similar  tran¬ 
sient  in  the  PE  signal  of  sodium  resorufin  in  dimethyl- 
sulfoxide.  Also  shown  in  Fig.  2(a)  is  the  best  fit  of  the 
Bloch  equations  to  our  simulated  results.  Clearly,  the 
Bloch  equations  do  not  offer  an  adequate  description  of 
the  initial  transient. 

Finally,  Fig.  2(a)  gives  the  results  of  mapping  the 
present  simulation  onto  the  Brownian  oscillator  model 
(6l.  These  results  closely  match  those  of  Eq.  (2).  The 
Brownian  oscillator  model  does  best  at  predicting  spec¬ 
troscopies  relating  to  ground-state  fluctuations,  such  as 
the  PE,  since  it  is  derived  by  an  expansion  about  the 
ground-state  dynamics.  Figure  2(b)  displays  lnS(f|). 
S(fi)  shows  a  pronounced  nonexponential  decay  at  long 
times,  in  contrast  with  the  exponential  decay  predicted  by 
the  Bloch  equations.  The  failure  of  the  Bloch  equations 
is  due  to  the  lack  of  time  scale  separation  in  C(r). 

Figure  3(a)  shows  the  calculated  3PE  signal  The 
<:“0  section  of  the  3PE  corresponds  to  the  PE.  The  3PE 
clearly  shows  the  vibrational  nature  of  the  solvent  dynam¬ 
ics  during  the  t ;  period,  in  contrast  with  the  PE.  Pump- 
probe  spectroscopy  (PP)  is  closely  related  to  the  3PE. 
We  consider  impulsive  PP,  where  the  pump  and  probe  are 
ultrashort  and  the  resulting  emission  is  spectrally  dis- 


FIG.  3.  (a)  The  three-pulse  photon  echo  signal,  (b)  The  im¬ 
pulsive  pump-probe  signal. 
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FIG.  4  The  nonlinear  polarization  |/>(,,(/i,/;.ri)|;.  The 
color  map  used  is  also  shown.  The  values  of  | /** "| '  were 
mapped  into  colors  with  the  color  map  shown. 


persed.  The  signal  is 

Spp(o}),t 2)  “Re  f  “0 )dt  ? . 

•J  0 

Figure  3(b)  shows  SVp  for  the  present  system.  At  /;“0, 
Spp  is  the  same  as  the  ordinary  absorption  line  shape.  As 
t>  increases,  the  excited-state  term  shifts  to  the  red.  lead¬ 
ing  to  a  broadening  and  a  growing  asymmetry  in  the 
spectrum.  The  Stokes  shift  of  the  excited-state  term  is 
not  predicted  by  the  Bloch  equations  or  by  the  stochastic 
model. 

PP  and  PE  are  complementary  special  cases  of  the 
3PE.  The  3PE  probes  the  full  nonlinear  polarization 
1T2T3),  whereas  the  PP  depends  only  upon  the 
G“0  plane  of  P(i\  and  the  PE  depends  onl>  on  the 
f2*“0  plane.  In  physical  terms,  the  PE  probes  the  decay 
and  rephasing  of  a  coherence,  whereas  PP  measures  the 
dynamics  of  an  excited-state  population.  The  3PE.  how¬ 
ever,  contains  both  the  coherence  decay  and  rephasing 
and  the  population  dynamics. 

Nonlinear  spectroscopies  measure  various  aspects  of 
P{})U\,ty*t })•  A  complete  measurement  of  Pty)  as  a 
function  of  its  three  time  arguments  could  be  done  by  us¬ 
ing  three  ultrashort  pulses  followed  by  a  fast  detection  of 
the  signal.  Figure  4  shows  the  complete  |/*<3>(r  2./ 1 ) 1 2. 


We  first  consider  the  /:  -0  plane  The  PE  signal  is  deter¬ 
mined  by  Piy)  in  this  plane.  P<}>  is  seen  to  show  an  echo 
at  times  / j“/i,  which  means  that  rephasing  of  the  initial 
coherence  occurs.  The  / 1  “0  plane  of  Pn>  corresponds  to 
the  PP  signal.  As  r:  is  increased  Ptu  decays  faster  with 
tc  This  corresponds  to  the  broadening  of  the  PP  signal 
as  the  Stokes  shift  begins,  since  1 1  is  the  variable  associat¬ 
ed  with  the  probe  absorption  line  shape. 

Finally,  we  consider  P{"  for  nonzero  /j  and  r?.  As  t: 
is  increased  from  0.  relaxation  processes  prevent  rephas¬ 
ing  from  occurring;  P'y>  is  less  peaked  about 
Only  at  very  short  times  can  the  thermal  distribution  of 
solvent  environments  be  considered  inhomogeneous;  later 
spectral  diffusion  during  destroys  the  echo.  We  also 
note  that  the  vibrational  dynamics  of  ('(/)  is  reflected  in 
the  oscillations  of  Pn>. 

in  conclusion,  the  essential  features  of  recent  experi¬ 
ments  in  polar  liquids,  such  as  the  non-Bloch  decay  of  the 
PE  signal  and  dynamics  in  the  intermediate  modulation 
regime,  were  captured  in  the  present  simulation  We 
were  able  to  clarify  the  relation  between  various  spectro¬ 
scopic  probes  and  the  underlying  molecular  motion.  In 
particular,  we  found  that  PP  and  3PE  experiments  are 
more  sensitive  to  vibrational  dynamics  than  PE  experi¬ 
ments.  Finally,  our  analysis  of  P,y>(/ :,r ,)  shows  that 
spectral  diffusion  during  the  /;  interval  very  rapidly  de¬ 
stroys  the  echo  character  of  the  polarization 
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